Quasi-Solitons in Dissipative Systems and Exactly Solvable Lattice Models 
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A system of first-order differential-difference equations witli time lag describes the formation of 
density waves, called as quasi-solitons for dissipative systems in this paper. For co-moving density 
waves, the system reduces to some exactly solvable lattice models. We construct a shock-wave 
solution as well as one-quasi-soliton solution, and argue that there are pseudo-conserved quantities 
which characterize the formation of the co-moving waves. The simplest non-trivial one is given to 
discuss the presence of a cascade phenomena in relaxation process toward the pattern formation. 
PACS numbers: 05.45.Yv, 05.45.-a, 04.20.Jb, 02.60.Cb 
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Recently three of the present authors have found ex- 
act solutions H] of the first-order differential-difference 
equations g] 



in{t + T)^V[AXn{t)] 

= ^ -\- rj tanh 



Axnjt) - p 
2A 



(1) 



where r is a time lag and Axn{t) — Xn~.i{t) — Xn(t) the 
forward difference on discrete indices. This system ex- 
hibits temporal-spatial density patterns of constituent el- 
ements, and can be used as a one-dimensional model for a 
variety of the co-operative phenomena such as congested 
flows of traffic, flows of granular particles, and possibly 
some chemical as well as ecological problems. The peri- 
odic solutionsn we have obtained precisely describe the 
steady profile density waves observed in numerical simu- 
lations. The density waves are locally stable and behave 
as soliton-like excitations. They are however neither nec- 
essarily "solitary" nor independent each other. It will 
be suitable to distinguish those appeared in dissipative 
systems from the solitons associated with integrable sys- 
tems. We call the new objects quasi-solitons (QS's). 

It is remarkable that the system (|l|) may exhibit a kind 
of partial integrability upon the spatial-temporal pattern 
formation. The idea has been suggested by three of us 
PI and also by the authors of ref. [H, who have written 



*An exact solution has been found also by Hasebe, 
Nakayama and Sugiyama W. ft is a special case of our 
solutions, obtained by a Landen transformation of elliptic 
functions. 



eqs.O) in a form which has a lattice analog, the system 
of Hirota's self-dual LC circuit equations. 

The purpose of this letter is to establish the connec- 
tion between the system (||) and relevant lattice models, 
constructing some exact solutions for an infinite (— cxd < 
n < cx)) system. Based on the reduction to the Hirota 
lattice [gl, we construct a shock- wave solution and one- 
QS solution. To get an insight into QS-solutions, we also 
transform eqs.(|l|) into the Hirota-Satsuma lattice model 
H , an extended version of the Hirota's model. One of the 
most important implications of the lattice model corre- 
spondence is that one may construct pseudo-conserved 
quantities in the dissipative system which quantitatively 
characterize relaxation processes toward the pattern for- 
mation. 

The lattice model correspondence applies only to the 
density waves which have the following properties: 
(i) a common phase variable w(t) = vt — kn\ 
(ii) a linear dispersion relation 



V = k/{2T). 



(2) 



The system (Q) describing propagation of the co-moving 
waves with these properties can be reduced to exactly 
solvable lattice models. Actually, the connection with 
the Hirota lattice is readily seen by defining 



Mr,{t) = tanh [(Aa;„(t) - p)/2A] , 



(3) 



and using (g) to replace the time shift t ^ t — t hy the 
subscript shift n -^ n + 1/2: w{t — t) — vt — k{n -f- 1/2). 
Making the total degrees of freedom twice of those of 
the original system (|l|), one obtains the Hirota lattice 
equations as reduced equations 



(1 - Ml) ( 



M„_i -Al,^i 



(4) 



where Tc — A/rj is a constant. As for the basic variables, 
we may take an ansatz 



Xn{i) = A In (pn{t) + Ct — nh, 



(5) 



where h and C = V{h) are constants. M is expressed by 

(f's as 

M„ = 'f"-^-'P" , ^ ^ exp[(;i - p)M]. (6) 

We give some solutions using these variables. 

(1) Shock-wave solution is given by 

ifn = l+ exp{-2w + k), 
M„ = tanh(/c/2) tanhw, (7) 

where the parameter should satisfy 

h = p + Ak, — A: = tanh(fc/2). (8) 

2r 

The solution (|^) yields a kink in the velocity variable 

i;„ = Av [tanh(w - fc/2) - 1] + C, (9) 

which gives an interpolation between a uniform flow with 
the velocity C — 2Ai' at n = — cxo and that with C at 
n = oo. Since the transition between two flows takes 
infinite time, the solution (Q) describes an asymptotic 
trajectory ||] in phase space. 

(2) One-QS solution takes of the form 



'Pn = 



1 + cx.p{-2w + k-2S) 
1 + exp{-2w + k + 2d)' 



M„ 



2C 



C + 1 (C + 1) 

Qg sinh(2(5) sinh fc 



,QS 

2 '^n ' 



cosh w + sinh d 



where 



k T AC 



sinhfc Tc (C + 1)^ ' 
tanh(/c/2) = -^—^ tanh(2(5). 

This solution describes a pair of kink and anti-kink 

in — Av [tanh(w — fc/2 + 5) — tanh(u) — fc/2 — 5) 
+C. 



(10) 

(11) 

(12) 



(13) 



We have obtained above a shock wave solution as well 
as one-QS solution. To get a further insight into QS so- 
lutions, we rewrite the reduced equations as those of a 
model of Hirota and Satsuma, who constructed exact so- 
lutions from those of the Toda lattice using the Backlund 
transformation technique. We show that this formalism 



gives the same one-QS solution as that of (|10D. To this 
end, we define new variables by 



2C 



u„ = Af„ - 



C-1 

C + 1' 



(C + 1)^ 

and a rescaled time variable 

2C ,_ 

Tcic+ir ■ 

These obey then the reduced equations 



ds 



1 : M 



C + 1 " (C+1)^ 



(14) 



(15) 



(16) 



The system ( |16D is exactly the same as the (u, v) system 
of the Hirota-Satsuma lattice. 



dun 
ds 

dVn 

ds 



a + PlUn + 



- 132V n 



{Pi - l)ul 



4a 



a(/3| - l)vl 



(w„_i -«„+!), 
(u„_i -w„+i)- 



with a = 1, Pi = 132 = -(C-1)/(C+1) and u„ = u„. The 
V variables are redundant here. In the Hirota-Satsuma 
construction, exact solutions are given by 



(17) 



Mn = -J- In ( /„//„ ) , 



where both /„ and /„ are solutions of the Toda lattice 
equations: 






1 + 312 11^/" = 



/n+l./n- 
J n 



(18) 



The variables u„ are suitable for solving the reduced 
equations subject to the boundary conditions Uqo = 
u_co = 0, which describe a uniform flow with Aa;„ = h 
at spatial infinities. As for one-soliton solution, we may 
take 



/„ = 1 + exp(2M/ + (/.), 
/„ = 1 + exp(2M/-H^), 



(19) 



where 



w ^ 


Cls — kn, 


n^ 


sinhfc, 


ex.])(j> - 


(1 + coshfc) — sinhfc, 

C + 1 


exp0 = 


(1 + coshfc) + sinhfc. 

C + 1 


This leads to 






nsmh[{4>-(t))/2] 



cosh2[W^ + {<j) + (/))/4)] + sinh2[(0 - 0)/4] ' 



(20) 
(21) 

(22) 



which should be compared with u'^^ in ( [lO| ) . Consistency 
in the dispersion relations requires ils — fct/(2T), which 
yields the condition ( [ll|). One obtains also 5 = {(j) — 
(/))/4, which leads to (|12|). Therefore, after making a 
constant shift of the time variable, two functions u'^^ 
and Un become exactly the same. 

It is well known that an exactly solvable lattice model 
has the same number of conserved quantities as those of 
the degrees of freedom of the model. The lattice cor- 
respondence discussed above provides us with pseudo- 
conserved quantities, which are generically time depen- 
dent, but become conserved quantities upon the forma- 
tion of spatial-temporal patterns of density waves. We 
may construct them from conserved quantities given by 
Wadati M for the Hirota lattice. For a finite system of p) 
subject to the periodic boundary condition Xi = Xi+N, 
the simplest one is given by a quadratic form of the M 
variables: 

N 

K{t) - Y. i^Mt) {Mn{t + r) -f Mn+i{t + t))] . (23) 

n=l 

It is not a conserved quantity, but after the co-moving 
density waves are generated, it reduces to 



N 



K,{t) = ^ [A'Ut) (Af„_i (t) + M„+i (i) 



(24) 



Using (Q), one can see that the reduced K actually con- 
serves, dIC/dt = 0. Thus K{t) characterizes relaxation 
processes toward exact periodic solutions which act as 
"attractors" of the system. In Fig.l, we show the time 
evolution of K{t) obtained from a numerical simulation 
of (11^) for iV = 20. The initial state is an almost uniform 
flow with K ^ 0.5, the first plateau in the figure. The 
uniform flow with small perturbation becomes unstable, 
and there appear density patterns with regions of high 
density where the constituent elements move slowly, and 
low-density regions where the velocities of elements are 
high. The high density regions of congested flows are 
viewed as clusters on the circuit. The final state with 
the maximum value for K r^ 12.6 is described by the ex- 
act periodic solution which contains one cluster. This is 
observed as the fourth plateau. In between, there appear 
the second and third plateaus of K ^ 7.7 and K ^ 10.1 
. These are the intermediate states well approximated 
by exact solutions with three and two clusters. There- 
fore, the relaxation process starting from a perturbed 
uniform flow exhibits a cascade decay via multi-cluster 
states. Our pseudo-conserved quantity K{t) quantita- 
tively characterizes this cascade phenomena. There exist 
other pseudo-conserved quantities which may take dif- 
ferent values for different pattern of flows, uniform or 
congested flows. Obviously we may construct inflnitely 
many kinds of such quantities with the same nature for 
an infinite system (— c» < n < cxd). 
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FIG. 1. The time evolution of the pseudo-conserved quan- 
tity K{t) in a simulation. 

We make several comments: 

(a) Our dissipative system is shown to be related to inte- 
grable systems and some solutions to the latter systems 
are translated in the former as congested flows. This 
connection, however, does not imply the stability of the 
solutions in the original dissipative system: they are re- 
lated only after patterns are formed. So the stability of 
the solutions should be studied in the original system. 
Actually we have confirmed that conditions in eq.(Q2) 
allow us to choose parameters to form typical stable pat- 
terns for congested flows observed in simulations M. A 
shock-wave solution may be understood as the extreme 
case of one-QS solution with (5 ^ oo, so the stability of 
the latter implies the same nature of the former. There- 
fore the solutions discussed in this paper are stable for 
suitable range of parameters. 

(b) In simulations we have observed the initial condition 
dependence; a flow might develop to a congested flow 
or a uniform flow depending on its initial condition. If 
we could label initial conditions with some of pseudo- 
conserved quantities, it would be a very interesting ap- 
plication of the quantities. These applications will be 
discussed elsewhere. 

(c) The dispersion relation ( H) was flrst given by 
Whitham for Newell's model ||l0| of traffic flows. It is 
a crucial condition which makes a bridge between the dis- 
sipative system (|l|) and exactly solvable lattice models. 
We showed that this relation is a kind of the intcgrability 
condition under which the system eq. (nh admits periodic 
wave solutions [^]. Furthermore, Hasebe, Nakayama and 
Sugiyama observed in numerical simulations that the re- 
lation in eq. (g) holds irrespective of details of the function 
V{Ax) B . It is therefore a mysterious but non-trivial dy- 
namical relation, whose derivation is challenging. 

(d) It is difficult to find a two-QS solution of the system 
(|l|). The co-moving limit of the Toda lattice solutions 
in which two-solitons propagate with a same velocity is 
not known. Therefore, one may not construct a two-QS 
solution from the Toda lattice solutions. It may not ex- 
clude, however, the existence of such solution. Actually, 
there exist periodic solutions of eq.(|^) for finite as well 



as infinite systems which do not solve the Toda equa- 
tions. To find a two-QS solution, one has directly to 
solve the reduced equations (|j) or (|l^) without referring 
to the Toda lattice solutions. Analysis of the co- moving 
limit in exact solutions obtained by Ablowitz and Ladik 
||ll| may give another insight into this problem. In this 
paper, we have restricted ourselves to exact solutions of 
co-moving waves. We do not know if multi-QS's with dif- 
ferent velocities, observed in numerical simulations, have 
corresponding analytic expressions. Discovery of such so- 
lutions, if any, is also challenging. 
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